N=10000;
ABC=[-1,sqrt(2)*1i,1] ;  %三个初始点坐标
Z=1i;
v=rand(N,1);    %生成0到1之间等概率分布的数
for kk=2:N
    if v(kk)<1/3
        Z(kk)=0.5*Z(kk-1)+0.5*ABC(1);
    elseif v(kk)<2/3
        Z(kk)=0.5*Z(kk-1)+0.5*ABC(2);
    else
        Z(kk)=0.5*Z(kk-1)+0.5*ABC(3);
    end
end
plot(Z,'.','MarkerSize',4)